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Abstract 
A computer program has been developed to expedite the calcu- 

lations outlined in the report entitled, "An Approximate Method 
for Selecting Optimum Drop Sizes" by R. Ekblad. The derivation 

of the equations upon which the program is based, together with 
the program listing and typical results, are presented. A second 
program entitled, "Minimum Drop Size Calculation", has been writ- 
ten to facilitate the calculation of the minimum lethal spray drop 


size. Comments on the relevant literature are also made. 
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Introduction 

The spraying of forest insects with an insecticide dispersed 
from an airplane or a helicopter is a very common method for 
controlling insect damage to tree stands. Several parameters 
are available for controlling the effectiveness of the spray. 
Some of these are: the distribution of the drop sizes, the wind 
velocity, the air temperature, the humidity, etc. These parame- 
ters are not independent in their action and the variation in the 
spray effectiveness due to a variation in these parameters is 
not readily predicted. In an attempt to analyze this problem, 

R. Ekblad has developed a method of analysis which is described 
in his paper entitled, "An Approximate Method for Selecting 
Optimum Drop Size", (ref. 1). It is assumed that the reader is 
familiar with the eee this report and consequently it 
will not be discussed in complete detail here. It is the purpose 
of this report to describe the development of a computer program 
based upon Ekblad's analysis. 

The objective of Ekblad's report is to describe a method for 
selecting the optimal spray drop diameter distribution given a 
characteristic target size, the physical characteristics of the 
air; such as the density, viscosity and velocity, together with 
the density and toxicity of the spray. The aoremeit She size 
depends upon the effectiveness of the spray which is estimated in 
terms of a base spray intensity. The base spray intensity is the 
number of drops of a given diameter appearing on a surface of one 
square a otaacer due to a total spray volume which is usually 
measured in macs per hectare. In this work, total spray volume 


will be measured in gallons per acre. 


A measure of the effectiveness of a spray is the probability 
of kill which is defined to be the total number of drops of a 
specified base spray intensity impacting a given target resulting 
inva "kill" of that target. ~Thus, by calculatingetiesprObabe ri. 
of kill for sprays with different distributions of drop sizes and 
the same base Aree intensitv, that is tne same total spray volume, 
it is possible to determine the optimum drop size for the specified 
base spray ares eee A target is said to be killed when it has 
received a minimum lethal dose of an insecticide. The minimum 
lethal dose is described in terms of a percent kill and is commonly 
denoted by LD, where x is the percent killed. Thus, LDg9 
Genotes a 90% kill. The lethal dose cepends upon the toxicity 
and concentration cf the insecticide. It is assumed that if a 
spray drop is of lethal size and touches a target, the target has 
been killed. The optimal drop size is determined by maximizing 
the total probability of kill as a function of drop size distribu- 
TON. 

The work that follows presents the analysis necessary to derive 
equations upon which the computer programs are based. Also included 
are the computer programs and their documentation. Two programs 
have been written. The first is entitled, "Optimum Spray Drop 
Calculation" and the second is entitled, "Minimum Drop Size Calcu- 
Lat LOn ar. 

The programs have been designed for ease of use by an analyst 
accessing the computer from a remote terminal in the time share mode. 


This programming style has been chosen because the program is to be 


pepntas 


an investigative tool. In addition, such a design makes the program 


easier to use in a field station since only a simple terminal is 


required as opposed to the necessity of having a remote job entry 


terminal if the batch mode of operation is used. 


Several computer runs have been included to illustrate the 


program 
Made in 


i036) 


program 


results. The selection of parameters used in the runs was 
accoLoawitchne the contract monitor. 
the following discussion frequent reference will be made to 


line numbers. This should enable an easier conclusion of 


the discussion with the program which is given in figures 9a-9}j. 


Figures 


VON. 


5, 6a-6c, 7a, 7b and 8 are flowcharts depicting the organiza- 


the program. 
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Technical Approach 


In order to determine the probability of kill for a spray con- 
taining a distribution of particle sizes, the distribution is char- 
acterized by a set of diameter classes and the number of particles 
of the spray in each class. The probability of kill is calculated 
for each diameter class for the same total spray intensity assuming 
that the spray consists only of drops of a constant diameter. The 
probability of kill for a distribution of drops is then calculated 
by taking a weighted sum of these probabilities of kill. The weights 
are the fraction of the total spray mass in each size category. Such 
a weighting is equivalent to weighting by fraction of volume of spray 
in each size class. 

The determination of the probability of kill for a specified 
diameter class requires the determination of the nose efficiency 
of a collection of particles all of which are of a constant diameter. 
The work of Brun et al (ref. 2) gives the impact efficiency of a 
collection of particles of equal size impinging upon a right 
circular cylinder. It is assumed the particles are distributed 
homogeneously throughout the air stream which is moving in a direc- 
tion perpendicular to the axis of the cylinder. The impact effi- 
ciency is the reciprocal of the ratio of the number of particles 
contained in the cross-sectional area swept out by the cylinder 
to the number of particles in the section which actually impinges 
upon the cylinder. Brun et al give the impact efficiency, EF, aS a 
function of two dimensionless parameters K and ¢, i.e. E=E(K, p ). 


K and p are defined as: 


x = 3.50824 x 107-12 Pw at u 


Al 


and 


p = 26.40 JES ee 


Af 


Koy is = density of air and water respectively in lbs/£t> 


a = radius of drop in microns 

U = velocity of air in miles per hour 

L = target radius in feet 

A4= air viscosity in lbs/(ft.sec.). 
Their results are presented in table 1 and appear in the program 
injplines/ tel 310eto> 1360 Mwher el x@andiyem denote ssw and p respectively. 

It is assumed that the target is a larva or pine needle whose 

shape may be approximated by a right circular cylinder of given 
diameter. From data on the toxicity of insecticides to forest 
insects, (ref. 3 and ref. 4), the desired dose may be obtained. 


The dose, or dosage, is the intensity of the application of the 
insecticide and is the ratio of the amount of applied pure insec- 
ticide to the body weight of the insect. The dose is usually 
expressed in units of micrograms of pure insecticide per gram of 
insect body weight or in units of ounces of pure insecticide 
applied per acre. The ratio of the volume of a lethal drop to the 


volume of the drop used in the impact efficiency calculation gives the kill ratio of 


the drop. The distribution of drops in the spray is characterized 
by the number of drops belonging to each diameter class of the 

spray. A diameter class is an interval resulting from the partition- 
ing of the total range of drop sizes in the spray. The product of 
the number of drops of a given diameter, the corresponding impact 
efficiency and the kill ratio gives the probability of kill for the 


number of drops of the given diameter in the diameter class. 


The sum of these probabilities for all the drop size classes in 
the spray is the probability of kill for the entire drop distri- 
bution. These calculations are accomplished in lines 4635 to 4680 
of the program. 

This method of estimating the total probability of kill 
omits many factors which are known to affect the determination 
of a quantitative estimation of the effectiveness of a spray. 
Some of the effects which are not accounted for in the analysis 
are the effect of evaporation, Of coagulation and) of £1] teringsave 
to the canopy or foliage of the tree. The role of turbulence or 
local air motion has also been ignored. The inclusion of these 
effects is not simple because both the experimental techniques 
and the theoretical analysis are formidable. However, it is possi- 
ble on the basis of what seems reasonable assumptions, to derive 
approximations for these effects. Future work will concentrate 
on developing such approximations and including them in the program. 

There are two common experimental techniques for the determi- 
nation of lethal doses. The first technique, called the topical 
application method, consists of individual carefully-measured appli- 
cations of the insecticide followed by a determination of the percent 
kill resulting from the application. The technique gives the dose 
necessary to kill a specified percentage of the larva. The specified 
percentage of killed larva is denoted by LD, and the required dose is 
expressed in units of micrograms of pure insecticide per gram of body 


weight of the insect. For each LD, there is a required dose. 


The other technique, called the spray application method, 
consists in depositing the insecticide on the insect larva with 
the aid of a sprayer. The drop and density distribution from the 
Spray is accurately controlled and the subsequent insect mortality 
determined. The required dose corresponding to a given LD, is given 
in terms of ounces of pure insecticide per acre corresponding to an 
applied volume of spray which is specified in gallons per acre. 

For both of these empirically determined dose statistics it 
is possible to derive a minimum drop diameter due to a specified 
spray concentration corresponding to a desired LD, - This deriva- 


tion is given in another section of this report. 
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Impact Efficiency Subroutine 


This subroutine is the heart of the program and is essentially 
a table look-up to determine the impact efficiency, E. Eisa 
function of the two dimensionless variables, KK and x F 

Since E is a function of two variables and is known only at 
specified points in the (K,@) space, it is necessary to use 
double interpolation (ref. 6) to estimate E whenever the point 
(K ,@ ) does not correspond to a grid point. The interpolation 
technique employed in this work is linear in both variables and 
uses a logarithmic interpolation scheme to estimate a pair of 
weighting parameters, A, and By: These parameters are then 


used to evaluate E from the expression 


E = (1-A,) (1-By)E, j+B) (1-A))Es 544A) (1-By) Ei 47 5*AgB Fuad 541° (2.1) 


Enecaquation 2.1, Ei 5 denotes E(K i, ®,) anc @iandm jemare 
the indices of the point specifying the quadrangle of the interpo- 

lation table containing (K,@). This;ee(i; ))  Characterizes the 

quadrangle whose pairs of opposite vertices are Givens Dy sn (1,)) ang 
fie lyeemancdm(it1)))enand (1, j+1).. The values of i and 5 


are determined in lines 8310-8335. The weighting parameters are 


given by 


= mate coma a (2.2) 
1 
InX,,, — 1X, 
and 
ine Ome nee eo 
eet Cee Pee) Sere eee 
oT (2.3) 


1n eel - In( B ,+1) : 


Equations (2.2) and (2.3) appear in lines 8410 and 8415 and 
equation (2.1) is accomplished in lines 8480-8490. The notation 
X = hK and Y= D has been used. Since D is zero, the value 
lis added to ? in equation 2.3 \toO insure the avoiding jot 
attempting "to take the logarithm of 0". This procedure results 
in at most a very small error, since P46 2 L007 

It may be the case that either MK or a) , or both and 
D , are outside of the range of the table. In this event special 
action must be taken. If either » or fe) are less than zero, 
lines 8230 and 8235, in conjunction with lines 8260 and 8270 
respectively, cause the program to stop. Similarly, if d 775070007 
lines 8245 and 8290 cause the program to stop. In all three 
cases the reason for stopping the program is printed out. It can 
also happen that YK > 320. In this event, the impact efficiency 
is set equal to one providing the value of p lies in the table 
range. See lines 8240 and 8277. An examination of the graphical 
representation of E(K go ), aS given in reference 2, shows that, 


at most, only a very small error is incurred by setting the impact 


efficiency equal to zero for large values of K . Whenever the 


above condition obtains, the computer prints out a statement like: 


"FOR SIZE CLASS NO. 6 X=466.326 AND X > X(1%), THUS E=1." 
(See for example, page 16.2 where this statement is printed out 1l 
times for size classes 6 through 16). The program variable xX 


denotes the dimensionless variable : 


For values of (X,@) in the range of the table, the deter- 
mination of E for small values of K and @_ requires special 
interpolation formulae. The next paragraph describes these inter- 
polation formulae and the region of the table in which they are 
used. 

Lines 8355 and 8360 set the impact efficiency equal to zero 
yon nO ite 0. ls e)<9)0.25, lines 8385 and 8470 calculate 


the impact efficiency from the expression 


E = 0.178958E, j texp(exp[5( 4X9 -0.1))) - 2.7128} . (2.4) 


This expression was obtained by noting the extreme exponential-like 
behavior of E with K whenever 0.1 < K < 0.25. If K > 0.25 
and p < 5, the variation of the impact efficiency with er) is 
negligible and hence the weighting factor Bl is set equal to 
zero. The calculation of E for this range of K and a) is 
given in lines 8440 to 8445 and line 8390 determines whether 


OUenot Dp hsein the required range. 


The flowchart of the entire subroutine is shown in figures 
favand sb. 

In order that the reader may more easily correlate the pre- 
vious development with the equations used in the program, the 


following table of notation is presented. 


Mathematical BASIC 
Notation Notation 
xn x 
Pp 5 
Ay Al 
By | Bl 
ak i 
3 J 
Bing ECL Ss) 


In terms of the BASIC programming notation, equation 2.1 


appears as: 


8480 LET El 


(1-A1) A (1-81) *E CL Oe teeR Cl Al) oer (sles) aie) 


8485 LET E2 


I 


AL* 1 —BIAB GL? ieee LSB 1 eed te) 


S490 P9 LETS Se= ie EZ 


Equations, 2/2.cand (2°53 are written as: 


8410 LET Al (LOG (X) - LOG(X(I)))/(LOG(X(I+1)) - LOG(X(I))) 


I 


9415. ©LETIBI COG Gia GOOG Gy, (J )rFL )))*)?7 (L0G (Y (J+1))" — LOG (Y (J) +1) ) 


and equation 2.4 in the program appears as 


Baye Bb Det. 1.6 950 eis ead) (EXP (EXP(5*(X—o1)'))) -° 2.7128). 


Calculation of the Probability of Kill Per Class 
Using the ASCAS Input Data 

The comparison of the effectiveness of different spray dis- 
tributions requires a determination of the drop distribution in a 
spray of prescribed intensity. The intensity of the spray is 
measured in units of liters per hectare or gallons per acre. An 
actual drop distribution is obtained from an experiment which 
consists in recording drop stain sizes appearing on sets of 
cards placed at prescribed locations in the spray rear The 
reading and the recording of the drop stain sizes is accomplished 
with the aid of an instrument and the ASCAS computer program 
converts the recorded data to actual drop size data. This sec- 
tion of the report describes the calculation of the probability 
of kill due to a spray of specified intensity and an experimen- 
tally determined drop stain size distribution furnished by the 
ASCAS computer program. 

The probability of kill for a given spray intensity is the 
sum of the probabilities of kill due each drop size class. The 
probability of kill corresponding to a size class is the number 
of drops in the size class, (N3) 5, times the impact efficiency 
corresponding to that drop size and the appropriate meteorological 
conditions during the spray. This assumes that every drop im- 
pinging upon the target is a lethal drop. If a minimum lethal 
drop size is required a (sees lines) 385) @thesp:Obebr tty moni beac 


obtained from the assumption that all drops are lethal must be 


reduced in proportion to the number of drops required to consti- 
ture a lethal drop. See lines 4635-4640. 

The total number of drops in the size class corresponding to 
a specified drop stain distribution will be denoted by (A7);. 
(A7), is determined by calculating the number of drops, (P7)., 
Aes FERS spray corresponding to the size class, and then weighting 
this number in accord with the fraction of the drop mass in the 
ee size class which fraction is obtained from the ASCAS drop 
stain data. (P7) ; is obtained by dividing the spray intensity 
by the volume of the average size drop in the drop size class. 
The diameter of the average drop is obtained from the ASCAS data 
in the following way. 

The ASCAS program gives the number of drops (N3) 5 ina 
given drop diameter class. The class boundaries are denoted by 
(S3); and are measured in microns. Thus, (N3) 5 denotes the 
number of drops in the size class whose inner drop stain diameter 
is (S3);_,+1 and whose outer drop stain diameter is (S3);. 

The size class boundaries are prescribed in terms of drop stain 

diameters rather than drop diameters since that is the way the 

boundaries are obtained from the automated reading instrument. 
th 


The diameter of the average drop stain in the i size class, 


(D3); is obtained from 


(55), 


i OID. 


($3). +1 
b, 





I. /p3)* = Z) 
Aa (DY, (Z { ($3), ~[ (53), +f 







Hence 
3 3 
(D 3) = ($3); - L(S3),, +L] 
L = . 
3{ (53), [sx ‘If 
where (D3) ; is measured in microns. In these expressions the 


term (S3),;_,+l is used rather than the term (S3); SIneCegiue 


eel 
counting mechanism reads in units of microns. Because the lower 
limit of integrationeis) 10, nuthesexpression for (D3), appears 


in the simple form 


(D3) = 


En 


(S3) 5 


The above calculations are accomplished in lines 3525 to 3550. 
The conversion of the average drop stain diameter to the average 
drop diameter, (R3)5, requires the use of the magnification 
factors appropriate to the particular set of ASCAS data. The 


conversion equation is 
3 = AT +e Bre (D3 )e-ct CF (D3)- 


where A, B, C and F are the magnification factors. Line 4545 


accomplishes this conversion. 


(N7) . is the number of drops of diameter (R3) 5 in an 
amount of spray covering one square centimeter resulting froma 


mono-dispersal spray of intensity Q liters per hectare. Thus, 


7), = £9 ¢ 200 


T (R3) 5 


)3 


or 
1.90761 x 10! 
(N7), = == PAX 
a 
(R3) 5 
where (R3) 5 is measured in microns. If the spray intensity is 


given as Q gallons per acre, the number of drops per square 
centimeter due to such a spray intensity, assuming a mono-disperse 
spray, is 


8 
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(R3) 3 


This expression is derived from the preceding expression by noting 
that one gallon per acre is equivalent to 9.3536969 liters per 
hectare. It is this latter expression which is used in the computer 
program because it is assumed that the spray intensity is expressed 


in gallons per acre. The expression appears in lines 4640 and 


4650 and the factor Q appears in line 4670. If the spray 
intensity is expressed in liters per hectare or some other set of 
units, the constants appearing in lines 4640 and 4650 would 
have to be changed accordingly. 

The fraction of the total mass of the spray contained in the 


see size class, (M4) 5, is calculated in the following way. 


Let (M3) 5 denote the mass of the average drop in the jen size 


class and let (M5) 5 denote the mass of all drops in that class. 
Then, 


(M5); = (N3); (M3); 


and the total mass in the spray, M4, is 


where S3 denotes the number of size classes. Hence, 
(M4) i = (M5) ;/M4, 


and the number of drops in the size class corresponding to a 


specified drop stain distribution is 


(A7) 5 = (N7) 5 (M4) 5 : 


iieitetseassumed that all drops are lethal, Pas the probability 
Ormekilil Of the drops in the sae class, is the number of drops 


impinging upon the target or 
Ps = E(A7) ; 


where E is the impact efficiency for that drop size. This 
equation appears in line 4670. The appearance of the intensity 
factor, Q, in line 4670 has been discussed above. 

To enable the reader to better understand the preceding work 
the following example is presented. The basic data is taken from 
an ASCAS print aut of experimental data obtained from the Beaverhead/ 
Gallatin Trial #3. This data corresponds to the data used in the 
attached program. The magnification factors are A=C=0, B=0.5556 
and F=l. The ASCAS program provides as input, the stain diameter 
class boundaries, (S3) 5, the mass of the average drop in the 
class, (M3) 5, and the number of drops in the class (N3);. This 
data is listed in tabular form in table #2 in columns l, 2, 5 and 6. 
Columns 3, 4, 7, 8 and 9 contain the results of the calculations 
corresponding to the respective equations in this section. The 
eatculation of (N7) 5 assumes that Q=l1. The tabular results 
should also be compared to the program results presented in 
Run amonapages: -16.13,to 716.16. In the program, the quantities 
labeled R3(Z) correspond to the quantities (R3) 5, the N7(Z) 
correspond to the (N7) 5 and the A7(Z) correspond to the (A7),. 


The output listed on pages 16.14 to 16.16, is an example of the 


auxiliary detailed output that may be obtained from the program 
if a 1 is typed in response to the query given by line 1120 

of the program. This ability to obtain detailed auxiliary 
results was provided to assist in the debugging and understanding 


of the program. 


Calculation of Number of Drops/Class From Cumulative 
Frequency Input Format 

This section of the report describes the method for obtaining 
the number of drops in each size class from the cumulative frequency 
vs. diameter curve. A size class consists of drops whose diameters 
lie in a given range. The range will be defined by (D5,D;,)) 
where D,; is the minimum diameter and Dial is the maximum diameter 
of the drop. Drops whose diameters lie in this range will be said 
to belong to the caae class. 


The cumulative frequency data will be given as a sequence of 


points (D;,C;) where C; is the cumulative drop frequency corres- 


ponding to the gee class size. Thus, the percent of drops in 
each class is 
Pee em — le oe 
See line 1725. 
For the ay class size, let Py, nh; and Ry denote the 


percent of drops, the number of drops and the average drop diameter 
in the class respectively. Also let N denote the total number of 
drops in a spray of Q gallons per acre and Vi denote the volume 


of a drop of diameter R,- Then, 


) Aa Q (2) 


or 
1 ty 
1 , n;R = 19 : (3) 


But 


Nis =_p.N (4) 


and hence, 


T 3 a 
6: Np, R; aos Q e (5) 


Thus, the total number of drops, N is given by 


Q 


Tr —_— 
> ) } 3 ; (6) 
6 PR; 


nt 


Ne 


Since Q is in units of gallons per acre and it is desired to know 
the total number of drops falling on 1 sq. cm., the above equation 


is modified to read 


_ 1.78642 x 10%Q (7) 


3 
ee oer 


af 


where R, is measured in microns. 


In the program Py is denoted by P3(I) and R; is denoted 
by D3(I). Lines 1750-1765 calculate the denominator of equation 
(2) and line 1770 corresponds to equation (2) with Q=l. In 


terms of the previous development it is possible to write equation 


(4) as 


8 
a; = 1.78642 x 10° 4 be 


3 
ns PR; 
and the probability of kill for the eae Size class is 


8 
P. = EQ D LOO 4 Zax 0 ; (9) 


a 
3 
2 PR 


The correspondence of this development with the program is accom- 
plished with the aid of the identifications: C3(I)=C,, D3(Z)=R,;, 
M4 (Z)=(M4) and P(Z)=P.. It is then seen that lines 4650 to 


4670 are equivalent to 


1.78642 x 10° 


R 
1 


Dea E(M4) . { 


i Oman (10) 


Equations (9) and (10) are equivalent. This can be shown by 


noting that 


3 
rEPSaE : (iT) 


From line 1800 we obtain 


pr? 
(M4) 5 = (12) 
M7 
or 
PR? 
(M4) 5 = 3 . (13) 
Deeb ARs 
ie 
Substituting equation (13) into equation (10) gives 
equation (9). Consequently, equations (9) and (10) are equivalent. 


The reader may wonder about the absence of such factors as 
the density, the total number of drops, and other constant factors 
in equations (10), (11), (12) and (13). These factors are not 
needed because the desired weighting factor (M4) 5 is a ratio of 
two masses and the missing factors would all cancel out. This may 
be seen by noting that the mass of the total spray in the Aen 


class is 


3 
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where S is the density and N is the total number of drops. 


The total spray mass is then 


"P wor p.R? 
6 1 
and the weighting factor is 
PRS 
Oe 7 : 
EePaR 


Thus, line 1800 does indeed give the weighting factor. Further- 
more, there is no necessity to include the factor Tp N/6 in 

the calculations in lines 1760 and 1800. The percent conversion 
factors will also cancel and therefore, they also are not included. 
Thus, M7 is not the actual mass but rather a relative mass. 

For the cumulative frequency vs. drop diameter input format 
the drops are assumed to be spread continuously over the diameter 
classes. Hence, the calculation of the diameter of the drop of 
average size ina size class, (D3) 5, is different than the calcu- 
lation used to calculate the average diameter when the drop data © 
is specified by ASCAS data. The fundamental equation is 

(S3); 
a te D“dD 


T 
6 bel, = eae 


T Pe 


6 
ioe 
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Thus, lines — 1655 and 1660 are different than lines 3535 and 
3540. 

As an example of the development described in this section, 
we use the data in the attached program. The necessary initial 
data is shown in columns 2 and 4 of table #3. The entries in 


the fifth column are obtained from 


8 
(N7), = 218642 x 0 


3 
Rj 


and the entires appearing in the last column are given by 


(A7), = (M4), (N7); 


The quantities (N7) 5 are labeled N7(Z) in the program and the 
entries (A7) 5 are denoted by A7(Z). Run 10a, pages 16.21 through 
16.2494 lists’ ‘a’ "detailed? printout and “teis seen that thesentcires 


in the last two columns of table #3 correspond to the printed results. 


sya 
Calculation of the Minimum LD. Drop Diameter, d 


The minimum LD, drop diameter depends upon the mass of the 
insecticide in the drop, the body weight of the insect, the concen- 
tration of the insecticide in the spray and the density of the pure 
insecticide. 

The units expressing the insecticide dose, D, required to 
obtain a specified LD, are given as KB gm/gm body weight or as 
oz./facre. Because of the two different units of measurement, the 
calculation of the minimum drop diameter, d, necessary to produce 
a given LD, is accordingly different. Diem istrparteot 
this section presents the derivation of an expression for d assuming 
that D is expressed in Pgn/gm body weight. d will be deter- 
mined in microns. 

Let: d = minimum drop diameter in microns 

D = Required dose in al gm/gm body weight, 
w = Insect body weight in mg, 


€ = Percent concentration of pure insecticide by 
volume of spray, 


and of = Density of pure insecticide in gm/cm? 


The calculation of the desired drop diameter requires a 
determination of the mass, M, of pure insecticide that must be 
contained in the drop. Equating this expression to the mass of 
insecticide actually contained in a single drop gives an equation 
for the determination of the drop diameter. 


The desired dose per insect is 


Bwaexbel OMe pam or Dwexel One ye 


of pure insecticide. The mass of pure insecticide contained ina 


spray droplet is 


Bio) fF x 10l? , 


where ,/dam1s “inemicrons: 


Thus, by equating these two expressions 


ee 3a 2c Ai =o 
= a? (ze) d x 10 = aDWexetO 
Owe 
43 = & Dw x 10° 
TO et , 
Hence, 
d = 57.5884 Dw 
Cr : 


If the dose is given in oz./acre, the determination of the 
mass of insecticide that must be contained in a single drop is 
accomplished by assuming that the spray is spread evenly over the 
area covering the insects. It is further assumed that the shape 


of the insect can be approximated by a right circular cylinder 


whose length is 4, cm and whose diameter is ¥, mm. THUS; 


the area, A, exposed by the insect to the spray is A £ax10-tem?. 


If D is the dose, M is given as 
» sie 
M = 7.0053 x 10 Jobe D gm. 


Equating this expression to the mass of insecticide contained in a 


droplet gives 


aie OES pee LS) 
7 0053" 10 4, 4,p == 4° (755) oe x 10 


One 


6 
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Hence, 





The program for these calculations is attached and the flowchart 


shown in figure 8. 
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Optimal Spray Drop Calculation Program Description and Comments 

The program accepts as input the physical characteristics of 
the air and the drops, together with the aropedtstribution otea 
particular spray, and calculates a PLroObablivteyporekiv igdistrabutdon 
and a total probability of kill. In addition, provision is made 
for considering a minimum lethal size drop. 

The overall organization of the program is displayed ie the 
flowchart depicted in figure 5. The flowchart is self explanatory. 
Figures 6a, 6b and 6c are a more detailed flowchart of the program. 
Tnesimpacteelficiency subroutane is®shownelin’ figures’: 7a:and 7h. 

The numbers appearing in the flowcharts refer to the appropriate 
line numbers in the program. 

The programs have been written in modular fashion in order to 
make them readily portable to another computer and also to permit easy 
alterations or additions. This has resulted in a program which 
appears lengthy; however, the actual ‘heart' of the program is rela- 
tively short. The generous use of remark statements adds to the 
lengthy appearance of the program and should enable the reader to 
eaSily follow the flow of the calculations. Such a method of pro- 
gramming optimizes the portability of the program to other computers. 
Portability of programs between computers is a serious and costly 
problem and every effort has been made to minimize possible diffi- 
culties which may arise in the transference of the program. The 
author deliberately chose simplicity over 'tightness' and '‘elegance' 


in programming style. 


Ore 


Because the experimentally determined drop distribution data 
may be furnished by the ASCAS computer program or it may be given 
as a cumulative frequency versus drop diameter curve, the program 
provides the user with the ability to select the data format 
appropriate to the problem of interest. The ability to give the 
user the choice of either type of data format requires that both 
types of data be stored simultaneously in the program. Since 
input data is stored in the computer with the aid of a data stack, 
and the operation of a data stack does not permit completely 
random access to the data in the stack, additional programming 
is required to furnish the user the option of choosing the type 
Of datasLormat: 

The data appears in the data stack in the order in which the 
data appears in the program and the data is taken from the stack 
to be used in the program in the order in which it appears in the 
stack. This order may not be the same as the order in which it 
is desired to use the data in the program. Consequently, if it 
is desired to use the data in an order which is different than. 
the order in which the data is stored in the data stack, it is 
necessary to first read in all of the data even though the data | 
will not be used immediately, or even at all, in the program. 
Thus, in this program, data from both types of formats is first : 
read in and then a decision is made as to which data sets will 
be used in the particular program run. This decision is made with 


the aid of a switch, F3. The flowchart appearing in figure 6b 
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may appear unusual because of the indirect method; however, the 


program does provide the necessary flexibility to the user. 
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Data Format Comments 


The Optimum Drop Size Calculation program will accommodate the two input 
formats most frequently used in presenting drop distribution data. Drop distribution 
data may be in the form of a cumulative frequency vs. drop diameter curve or a 
set of points defining such a curve. Acceptance of this form of 
drop distribution data is accomplished in lines 1632-1636 and in 
lines 1702-1706. The other form of drop distribution data is 
that given by the ASCAS program. The program provides for this 


format in lines 3270-3274, 3340-3346, and 3440-3446. 


For the cumulative frequency vs. diameter format of data input, 
the actual droplet diameters are given as upper boundaries as opposed 


to the ASCAS format which gives stain diameter boundaries. The set 


of actual drop class limit diameters is read in (see line 1625) 

using the labeling variable S3(I). This is the same variable used 

to denote the stain diameter boundaries when reading from an ASCAS 
type format. Because the program uses the magnification factor 
formula, line 4545, to transform average stain diameters to average 
drop diameters, the values of the magnification factors, A, B, C 

and F must be suitably set whenever the S3(I) are actual droplet 
diameters (i.e. when using the cumulative frequency vs. drop diameter 
input format). This is done by setting A=C=0 and B=F=1; thus enabling 
line 4545 to merely relabel the S3(1I) to R3(1I). This modifica- 
tion of the magnification factors takes place whenever the cumulative 
frequency vs. drop diameter input format is used and permits both 


input formats to be accommodated in the same program. 


For the attached set of computer runs, the ASCAS data is 
taken from spray data obtained in the Beaverhead/Gallatin Trial #3 
experiment. The spray density was 1.011 gm/ml and the magnifica- 
tion factors were A=0, B=0.5556, C=0 and F=l. The cumulative 
frequency data is taken from a sample set of data provided by the 
contract technical monitor. Both of these sets of data are used 
for illustrative purposes only. If other drop distributions are to 
be used as a basis for the calculations, they must be entered into 
the data statements in accord with the preceding instructions. If 
the acquisition of such data were automated, the modular style of 
the programming would easily permit the DATA input section of the 


program to be altered to accommodate the automation. 


The User Guide section explains in greater detail the DATA 
input. 

The Minimum Drop Size Calculation program accepts as input the 
physical and entomological parameters required to determine the min- 
imum lethal drop diameter. These parameters are listed in the sec- 
tion describing the derivation of the minimum drop size and are 
called for, in proper order, by the program as it requires them. 

The notation at the beginning of the program listing gives the units 
of the parameters and the desired drop diameter is calculated in 
microns. A flowchart of the program appears in figure 8 and the 


operation of the program is provided by the program as it is run. 
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USER GUIDE 


Data Input 

Before the user runs the program, it is necessary to enter as 
DATA the detailed description of the particular spray. As stated in 
section 7, there are two formats used in the detailed description of the 
SOLay.ee ne Lirst format. is that which is produced by the ASCAS 
computer program. From that program the stain class size diameters 
in microns, S3(1I), may be obtained and are entered as DATA in 
ness 7 Oto 3274. The ASCAS program also gives the number of drops in each size 
Class which are entered as DATA in lines 3340 to 3346 together 
with the mass of a drop of average size in each class which is 
entered as DATA in lines 3440 to 3446. The number of classes, or 
SzemealecOCles sao ,.0.1S entered by the operator as.input in. line 
1597. The second type of format used in the description of a spray 
is a set of points obtained from a cumulative frequency vs. drop 
diameter curve. The drop diameter class boundaries are entered as 
DATA in lines 1632 to 1636 and the cumulative frequencies corres- 
ponding to each class boundary are entered as DATA in lines 1702 
to 1706. In the event that more than 20 size categories are required 
to describe the spray distribution, the dimension statements must 
be altered to accommodate the increase. Presently, up to 20 size 


classes can be accommodated. 


During the course of running the program, the following 
additional data is required. The density of the air, D2, in 
ipew/tee the viscosity ot the ain, 9 Vl, sinelbs./ tt aseesns 
the density of the spray droplet, Dl, in lbs. /ft. >; the diameter 
of the target, Cl, in inches; the velocity of thegair, S.Uljwaa 
miles per hour; and the base spray intensity, Q, in gallons per 
acre, must all be supplied by the operator. These six pieces of 
data are entered in lines 4270 and 4340 of the program and are 
called for by the program via printed instructions on the terminal 
to the operator. 

If the ASCAS input format epbhen fie selected, the magnifica- 
tion factors A, B, C and F must be supplied by the operator. They 
are called for by lines 3155 and 3165. Since the ASCAS computer 
printout provides the drop diameter and the percent mass by class, 
the total mass, and the number of drops per size class, it would 
have been possible to directly enter such data and, as a result save 
some calculational effort. However, it may be the case that only 
raw stain data and the attendant magnification factors are available, 
and hence provision was made for using such data. In addition, the 
extra calculation required because of the acceptance of this form 
of the data is very very small and consequently it was not thought 
worthwhile to omit the capability. 

There are several operator options available. The first option 
appears in line 1120 and allows the operator to obtain a detailed 
printout of the calculation of the essential quantities necessary 


to calculate the probability of kill for each size class. These 


calculations are performed in lines 4635 to 4672. The quantities 
printed out are listed in lines 4684 to 4687 and may be identi- 
fied with the aid of the notation listed at the beginning of the pro- 
gram. Examples of such detailed printouts are given in runs #7a and lO0a 
in the attached set of program results. The second option, which 
appears in line 1385, permits the specification of a minimum 
lethal drop diameter. This diameter must be expressed in microns. 
The third option permits the designation of either the ASCAS 
input format or the cumulative frequency vs. the drop diameter input 
format. The option appears in line 1580. The last option, line 4760, provides 
the user with the ability to run the program again, using all input 
variables and data as before, with the exception of the target 
diameter, the air velocity and the base spray intensity. These 
latter three parameters are varied most frequently and this option 
enables a parameter study involving the effects of these three 
variables to be readily carried out. The option is invoked by 
typing a 1 in response to the request indicated by line 4760. The 
output necessary to construct figures 3 and 4 of the Ekblad report 


is suppressed when this option is exercised. This output is not 


necessary Since the first running of the program provides the 
output. 

Since any combination of choices of these four options is 
available, the program is quite flexible and, it is hoped, will 
thus readily accommodate the user. Examples of the use of these 


options are presented in the attached set of program results. 


Data Output 


The results of the program are presented in the form of tables 
containing indices and headings for identification. Some typical 


results are shown in the attached set of program runs. 


Output of Results 
The program presents the results as a set of tables whose 
entries correspond to that displayed graphically in figures 3, 4, 6 


and 7 of the report by Ekblad. 


Program Results 

Results from several runs are given in order to illustrate 
the different output formats and to also indicate the potential of 
the program as an aid in studying the variation of the probability 
of kill with the parameters. For the Optimal Drop Size Program 
an index of program runs and resultsis provided on page 16 to enable the iden- 
tification of the values of the input parameters with the corresponding runs. 

The setting of the minimum lethal drop size, Ll, equal to 
zero corresponds to the assumption that every drop is a lethal drop. 
piemseeting OL Li-134.5 "and “Li=25/7.2 “corresponds to a topical 
application of mexacarbate to the ee and Bae instars respectively at the gp 
level. The choices L1=124.8 and 285.8 correspond to spray appli- 
cations of mexacarbate to the a and act instars respectively at the IDgp 
level. These insecticide data, together with the target diameter 
data, are for the spruce budworm and were obtained from the contract 
monitor. References 3 and 4 illustrate the format in which such 
data is presented. 

Runs 7a and 10a have the same set of parameters as runs 
7 and 10 respectively and the purpose of runs 7a and 10a is 
to illustrate the output produced from lines 4684-4687. This option 
was provided to enable easier checking of the calculation. 

On several of the runs the calculated value of the dimension- 
less parameter YW exceeded the allowed range. This indicated 
that the range of the table had been exceeded and consequently, the 
impact efficiency was set equal to unity. For those runs for which 
this occurred, the size class together with the value of K , 


(ca bleamexmeinethesprogram), was printed out. The exceeding of the 


range of the table usually occurs when the target radius is small 
and the air velocity is large. This is due to the fact that 
varies inversely as the target radius and directly as the air 
velocity. An illustration of the results furnished by the program 
when x has exceeded the allowable range is provided by Run l, 
page 16.2. The lower portion of the printout indicates that the 
range of K , for this particular run, had been exceeded for all 
size classes greater than the fifth size class. 

The figures referred to in the program listing and in the pro- 
gram output refer to the corresponding figures in the Lepontaay 
Ekblad. 

Results from the Minimum Drop Size Calculation program were 
selected to be illustrative of insecticide test data on the spruce 
budworms. The program is listed on page 15.11 ana typical sets of 
results are shown on Bag ese lon) 2eandml 5913) 

Figures 1, 2, 3 and 4 are graphical representations of a por- 
tion of each of the results obtained from runs 12-15 respectively 
of the Optimum Spray Drop Calculation Program. Since the area under the curve is 
a crude indicator of the overall effectiveness of the spray, it is quite evident that 
an increase in target diameter considerably reduces the effectiveness of the 
Spray. The graphs also show that there is relatively little 
change in the Spray effectiveness if the drop velocity is increased 
from one to six miles per hour. Certainly, for much higher veloci- 


ties it is not expected that this would be the case. 
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Literature Review 

A review and analysis of the relevant literature was made to 
enable the relating of the present work to that of others. In 
particular the work concerning the calculation of impact efficiency 
was reviewed. The original work of Langmuir, et al, ref. 6, used 
a mechanical differential analyzer to calculate the particle tra- 
jectories. Brun, et al, using the theory of Langmuir's as a foun- 
dation, obtained more accurate results by using a more accurate dif- 
ferential analyzer. They also showed that the effects of compressi- 
bility, for the air speeds of interest, were indeed negligible. 

Johnstone, et al, ref. 7, discussed the dispersion and deposi- 
tion of aerosols. However, their work was devoted primarily to small 
particles, i.e. those whose diameters were less than 100 microns. 
They used a diffusion model whose validity becomes more questionable 
as the particle size increases. In a review article entitled, 
"Filtration of Aerosols by Fibrous Media", Chen, ref. 8, summarizes 
the various parameters affecting the calculation of the impact effi- 
ciency. Some of the parameters are: inertia, Brownian motion (for 
very small particles), settling, diffusion and the effect of a collec- 
PvOnusOtepacticlessonsthe overall flow pattern. This article is a 
good review article up to 1955. 

Dorman wei naascnaptersentitied, Filtration”, sin *thestext edited 
by Davies, ref. 9, discusses the parameters affecting the calcula- 
tion of the impact efficiency. In addition to the many parameters 
considered by Chen, Dorman also considers sedimentation, electrical 


effects, packing density and clogging effects. Both Chen and 


Uy 7 


Dorman's articles have quite complete bibliographies. In the same 


| 


text, the chapter by Pich entitled, "Theory of Aerosol Filtration 

by Fibrous Membrane Filters" is very inclusive. The article contains 
a good summary of the comparison of experimental and theoretical 
results together with a very extensive bibliography. Davies, in 

the same text, considers the deposition from moving aerosols. His 


work contains only a small section on calculation methods for obtain- 





ing the impact efficiency. “The work’ of "Green and” Lane; ret wip, 
covers much of the same material, only in a more general 


setting. 














Many other articles were also surveyed. However, they were 
not found to be sufficiently closely related to the present work 
and so are not included in this discussion. Some other articles 
that may be of interest are listed as references 11, 12 and 13 


respectively. 


2. 


10- 


dj bales 


Le. 


The ts 


BEML a @h 


REFERENCES 


Ekblad, R., “An Approximate Method for Selecting Optimum 
Drop Size". Unpublished Report. U. S. Department of 
Agriculture Forest Service, Equipment Development Center, 
Missoula, Montana, 1976. 


Brun, R. J., Lewis, W., Perkins, P. J. and Serafini, J. S., 
"Impingement of Cloud Droplets on a Cylinder and Procedure 
for Measuring Liquid-Water Content and Droplet Sizes in 
Super-Cooled Clouds by Rotating Cylinder Method". NACA 

ENG oe 9.55, 


DON pe eee rOWl aS ead sand sRichmonda, C. E., “Insecticides 
Tested on New Tussock Moth". Bulletin #161, Pacific Southwest 
Range & Experiment Station, Forest Service, U. S. Department 
Of Agriculture. 1967. 


ROvemE SON em ci lete Nea ly.) LOOK, M., Lucas, B. A. and 
Lyon, R. L., "Toxicity of Selected Insecticides Applied to 
Western Spruce Budworm". Journal of Economic Entomology, 
Voie ONO: 0976. 


Carnahan, B., Luther, H. A. and Wilkes, J. D., "Applied 


Numerical Methods". J. Wiley, N. Y. 1969. 
Langmuir, I. and Blodgett, K. B., "A Mathematical Investigation 
of Water Droplet Trajectories". NACA TN 5418, Feb. 19, 1946. 


HOnNnsStOne, s.r.) WINnSChe, We &. and Smith, L. W., "The 
Dispersion and Deposition of Aerosols". Chemical Reviews, 
#445 April 1949. 


Ciena tee bt Ltrati1On ot Aerosols by Fibrous Media". 
Chemical Reviews, Vol. 55, June 1955. 


Davies, C. N., "Aerosol Science", Academic Press, N.Y. 1966. 


Green, H. L. and Lane, W. R., “Particulate Clouds, Dusts, 
Smokes and Mists", 2nd Edition. University Press, Belfast, 
Ireland, 1964. 


Golovin, M. N. and Putnam, A. A., “Inertial Impaction of Single 
Plenents =e fnd. Eng. Chem. Vol. 1, #4, 1962. 


Siti at eanawuoOrdan, M. L., “Mathematical and Graphical 
Interpretation of the Log-Normal Law for Particle Size Dis- 
trabution Analysis". Journal of Colloid Science, 19, 1964. 


Brooks, F. A., "The Drifting of Poisonous Dusts Applied by 
Airplanes and Land Rigs". Agric. Eng., Vol. 28, No. 6, June, 
od 





ifs Uh 


Summary 

The present work describes the rationale and the development of a computer 
program which implements the predictive technique of R. Ekblad for prescribing 
Optimal air dispersal of sprays. As stated in the technical approach section, several 
important factors have been omitted from the analysis and hence 
from the program. Both the author and Mr. Ekblad plan to extend 
the analysis to include these effects. 

An examination of the results of the program runs illustrates 
the flexibility of the analysis and the program. Since the intent 
of this work was to develop the computer program and to document 
this development, as well as to demonstrate the program results, 
no parameter study was attempted. Consequently, an analysis of the 
effect on the probability of kill due to a variation of the spray 


parameters will not be presented. 
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Sample Drop Distribution Calculation 
from Cumulative Frequency Data Input 
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Run No. 13 
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Total Prob, of Kill Index = 7.62 
Air Velocity = 6 mph 
Target Diameter = .03937 inches 


Probability of Kill Index 
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Run No. 14 
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Total Prob. of Kill Index = 4.55 
Air Velocity = 1 mph 
Target Diameter = .11811 inches 
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DETAILED FLOWCHART (continued) @) 
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DETAILED FLOWCHART (continued) 
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Impact Efficiency Subroutine (continued ) 
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Minimum Drop Size Calculation 
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REM 
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1610 REM LINES 1620-1626 READ IN DROP SIZE BOUNDARIES ¢C. F. 
1644 REM 

4614 REM 

14615 REM LROF SIZE BOUNDARIES IN LINES 1622-1626 MUST BE 
1616 REM CHANGED FOR EACH NEW SET CF EXPERIMENTAL DATA 
1617 REM 

4626 FOR I=1 TO Sz 

14625 READ SECT 
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4634 DATA 210.4, 250. 2, 284, 7, B42, 4, E97. 2, 458. 1 

14626 DATA Sid. 5.554. 7. 624. & 
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=330 NEXT I 

=340@ DATA SS&440e, 261277, 601465, 558674 
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Index of Optimal Drop Size Program Runs and Results 
For all runs the following variables were held constant: 
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18 Cum. Freq. 134.5 OS 957, 1 1 
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